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ABSTRACT 

We present an update to seven stars with long-period planets or planetary candidates using new 
and archival radial velocities from Keck-HIRES and literature velocities from other telescopes. Our 
updated analysis better constrains orbital parameters for these planets, four of which are known multi¬ 
planet systems. HD 24040 b and HD 183263 c are super-Jupiters with circular orbits and periods 
longer than 8 yr. We present a previously unseen linear trend in the residuals of HD 66428 indicative 
on an additional planetary companion. We confirm that GJ 849 is a multi-planet system and hnd a 
good orbital solution for the c component: it is a IMjup planet in a 15 yr orbit (the longest known for 
a planet orbiting an M dwarf). We update the HD 74156 double-planet system. We also announce 
the detection of HD 145934 6, a 2Mjup planet in a 7.5 yr orbit around a giant star. Two of our stars, 

HD 187123 and HD 217107, at present host the only known examples of systems comprising a hot 
Jupiter and a planet with a well constrained period > 5 yr, and with no evidence of giant planets in 
between. Our enlargement and improvement of long-period planet parameters will aid future analysis 
of origins, diversity, and evolution of planetary systems. 

Subject headings: planetary systems — stars: individual (HD 145934, GJ 849) — techniques: radial 
velocity 


1. INTRODUCTION 
1.1. Background 


The possibility of “Earth 2.0”, and especially another 
planet that hosts life, drives much of the search for exo¬ 
planets. As of 2014 November, efforts over two decades 
have uncovered more than 1400 planets and almost 4000 
planetary candidates (Han et al.|2014 Burke et al.|2014 


; exoplanets.org). The variety of discoveries, trorn lone 
Jupiter-mass planets in few-day orbits to packed systems 
with multiple planets that fit within Mercury’s orbit, 
raises a significant question as to the nature of our solar 
system: Are we unique? 

To search for analogs of the solar system, we target 
multi-planet systems and long-period giant planets, rem¬ 
iniscent of our own outer solar system. Because we seek 
planets with orbits of at least a few hundred days, the 
radial velocity ( RV) method of exoplanet detection is ad¬ 
vantageous (e.g.lWright fc Gaudi 2013 Butler et al. 1996 


Mayor fc Quelo^I995[). The RV method is the longest 


running, with rnuftiple surveys studying thousands of 
stars. In our study, we utilize up-to-date velocities from 
Keck observatory’s High Resolution Echelle Spectrome- 
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ter (HIRES; Vogt et al. 


19941 and complementary, pub¬ 


lished velocities from other telescopes, where available. 

So far, the number of planets discovered by RVs with 
periods greater than 1000 days is 103, only 16 of which 
have periods longer than 3000 days. This is in contrast 
to the 336 such planets with shorter periods]^ Our intent 
is to enlarge this sample of long-period planets to search 
for planetary systems with Jupiter analogs. Of these 103 
planets, 31 are in multi-planet systems. The study of 
multi-planet systems addresses planetary formation, mi¬ 
gration, and dynamics. Having a large sample can also 
contribute to the understanding of the evolution and life¬ 
time of stable planetary systems. Studies can examine 
the orbital eccentricities and perform dynamic simula¬ 
tions and probe migration. _ 

F or the purposes of this discussion, we follow [Wang et| 
[ah] (2012) and dehne a Jupiter analog as a planet with 
D > 8 yr, 4 > M sin(i) > O.bMj^p, and e < 0.3, but we 
also adopt an upper period limit of P < 16 yr. Of the 
conhr med RV planets , only 13 planets fit the above cri¬ 
teria ( Han et al.|2014 ). Another motivation for studying 
systerns with loiig-period Jupiter analogs is the role such 
a Jovian planet may play in the habitability of an Earth- 
like planet in the same system. Wetherill (19941 argued 
that Jupiter acts as a shield that deflects comets originat¬ 
ing from the Oort Gloud or Kuiper Be lt, protect i ng the 
inner solar system. Without Jupiter, Wetherill (19941 
suggested an increase in the frequency of cometary im¬ 
pacts on Earth by 1000 - 10,000 times the present-day 
value. Multi-planet systems serve not only as examples 
of planet-planet interaction but also as models for plane- 
comet dynamics. 


We follow [Han et al.j |2014| , who adopt a upper limit on min¬ 
im mass of 24 times the mass of Jupiter 
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1 . 2 . Plan 

Section gives an outline of the steps taken for char¬ 
acterizing the planetary systems. In Section we de¬ 
scribe the planetary systems orbiting seven stars. Each 
of these systems already has at least one planet known 
and exhibits RV residuals indicative of an outer com¬ 
panion. Additionally, each can have its planetary orbits 
signihcantly rehned with our new velocities from Keck, 
and in some cases we show that an outer, decade-long 
planetary orbit has finally completed. We present a sum¬ 
mary of the radial velocity data, mean uncertainties, and 
telescope offsets in Table Table lists the stellar pa¬ 
rameters of the target stars. Table lists the orbital 
parameters of the planets presented in this paper. We 
also present figures showing the RV curves and residuals 
for each system. Section presents an analysis of a new 
planet, HD 145934 b. We aiscuss our findings and future 
prospects in Section 


2. METHODOLOGY 

2 . 1 . Radial Velocity Sources and Analysis 

We combine previously published data from other tele¬ 
scopes to complement the time span and quantity of 


Planet Survey (iHoward et al. 112010 

lohnson et al . 112010 

Wright et al.|2Ull 1 tor many purposes, including as part 

of the ?70 survey (Howard et af. 

2009 

2010 

2011 a|b 

). At 


(Mayor fc Queloz||1995 ), several surveys were underway 


and actively monitoring stars for the signs o f planets (e.g. 


Cochran & Hatzes 1994 Fischer et alTl2014|). The discov- 


ery team for 51 Fega si used the ELUDIE spectrograph 
(Baranne et al. 1996), which was part of the Northern 
Extrasofar Ffanet Search until t he SOPHIE spectrograph 
(Bouchy & Sophie Team 2006) replaced it in 2006. The 
(JURALIE spectrograph (e.g. Queloz et al. 2000) was 
situated in Chile as part of the Southern Sky extrasolar 
Planet search Progra mme. It has been joined by HARPS 
(Mayor et al. 2003), also located in Chile. We make 
use of literature data from all four spectrographs in this 
work. Other data come from High Resoluti on Spectro - 
graph (HRS) of the Hobby-Eberly Telescope (|Tull|1998|), 
the Tull Spectrog raph at the 2.7-m telescope of McDon- 
ald Observatory (Tull et al. 1995 ), and the Hamilton 
spectrograph at Lick Observatory ( Vogt||1987 ). _ 

To analyze and fit the data, we use the [Wright fc| 
Howard (2009) RVLIN package written in IDL that natu¬ 
rally handies multiplanet systems using data from mul¬ 
tiple telescopes in systems where planet-planet interac¬ 
tions are negligible given the precision and the span of the 
observations. In this package, RV curves are described by 
both non-linear and linear parameters, and the package 
performs least-squares htting on them separately. The 
package uses a simple linear least-squares solution for 
the linear parameters, and the Levenberg-Marquardt al¬ 
gorithm for the nonlinear parameters. RVLIN supplies 
a sum-of-Keplerians model (plus optional secular trend 
and offsets between instruments) to MPFIT , the IDL im- 
plementation of the LM method developed by|Markwardt| 
f200^. 

We ht for these new planetary system as follows: (1) 
we collect Keck RVs before and after the 2004 HIRES 
upgrade separately to account for any (small) offsets be- 


tween the p re- and post-upgrade time series (e.g., Kane 
et al. 2014[); (2) we collect literature RVs for the sys- 


tem from other telescopes, if available (3) we use the 
published orbital parameter s (which we collec t from the 
Exoplanet Orbit Database; Han et al. 2014) as initial 
guesses for the planets’ orbits; (3) we use RVLIN to fit 
the system anew (with additional planets contributing 
five model parameters each, if necessary); (4) we use the 
reduced X^iXv) to describe the goodness of fit. 

We calculate most orbital para meter uncertainties us¬ 
ing BDOTTRAN ( Wang et al.||2012 ), which uses RVLIN and 
a bootstrapping method to compute the distribution of 
parameters consistent with the data. Because uncertain¬ 
ties can be highly non-Gaussian for planets with incom¬ 
plete orbits, we also examine the minimum surface in 
minimum mass-period space (Section^. 

For our hts, we choose values for the jitter (Wright 


2005 and references therein) that yield Xu values close 


to 1 ; usually we pick a value similar to the rms of the 
initial fit which does not incorporate jitter. If a star has 
data from several (more than three) instruments taken 
by multiple teams, we apply jitter on an instrument- 
by-instrument basis. To do so, we ran the fit with no 
assumed jitter, calculated for each instrument the stan¬ 
dard deviation of the residuals, and added that value in 
quadrature to the velocities. After that, we rerun the 
fit and that yielded the best-ht parameters. We utilized 
an instrument-by-instrument jitter for HD 24040, HD 
74156, and HD 217107. In general, we are confident in 
the relative instrumental uncertainties in the pre- and 
post-upgrade HIRES data, and we use a common jitter 
value for bothj^ Table lists the stellar parameters of 
host stars, and Tablelists the orbital parameters of the 
planets we discussed below. Corresponding RV plots and 
additional figures follow the text. 

Sun-like stars are known to have magnetic cycles with 
periods comparable to the period of Jupiter (Baliunas et| 


al.|1995|). A persistent concern in the hunt for long-term" 
W 


RV signals from Jupiter analogs has been that they might 
be mimicked by the effects of such magnetic cycles, which 
could alter convective patterns such that the magnitude 


vary with the stellar cycle dDravins 

1985 Walker et al. 

1995 |Deming et al.||1987 

Santos et al.||2UlU). 


A common way to check that magnetic effects are not 
responsible for RV variations is to measure correlations 
between the RVs and a ctivity indices such a s Ca II H 


fc K. Pre vious work by W right e t al. (2008), Santos et 
al. ( 2010| , and Lovis et af.ld'ioill) hnd that the observed 


activity-cycle-induced RV amplitudes are typically quite 
small (a few m s“^ or less), although there are sugges¬ 
tions that a few stars may show abnormally high levels 
of correlation (at the level of 10-20 m s“^). 

We have checked for activity cycles in these stars to 
see if they have similar periods and phases to the RV 
measurements. To do this, we have used the Ca II H 


^ The effect of jitter on the best fit values of an orbital solution 
is to give more even weight to points with different measurement 
uncertainties; in the cases of the well-detected planets we discuss 
in this work, the exact value of the jitter has very little effect on 
these best-fit values. Because we determine most of our parameter 
uncertainties via bootstrapping, our uncertainties are not strongly 
affected by our choice of jitter, and so there is no need to find the 
precise jitter value that yields x? = l-O- 
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fc K chromo spheric activity measurements from |Wright| 
et al. ( 2004| hereafter W04), Isaacson & Fischer (2010 


hereafter Ih'lO), and more recent measurements made us- 
ing the same data stream and pipeline as the latter work. 
For some stars, there appear to be calibration differences 
between the measurements published by W04 and those 
made using the IF 10 pipeline, necessitating a rescaling 
or application of an offset to one of the streams. This 
is most apparent in “flat activity” stars which show no 
variation but occasionally exhibit a large jump in activity 
level between the two data streams. 

In six of our stars, there is no appreci able activity vari - 
ation (i.e., they are “flat activity” stars, iSaar et al. 1998), 

m 


making it very unlikely that the large RV variations we 
see are due to solar-type activity cycles. The seventh 
star, HD 183263 does show a significant cycle, however. 
The W04 activity levels decrease from 2002 to 2004, and 
the IFIO show a continued decrease starting in late 2004, 
which bottoms out in a minimum in 2012. The actual 
velocities show a minimum in 2005 and a maximum in 
2012 , thus exhibiting a shorter period than the actual 
activity cycle. The negative correlation between 2005- 
2012 and the very high amplitude of the RV signals are 
inconsistent with typical s tars w ith RV-activity c orrela- 


tion seen m 


Wright et al. (2008 1 and described by Lovis 


et al. (2011 r it is thus very unlikely that any of the 
long-period signals we describe in this work are due to 
stellar magnetic activity cycles. 


2 . 2 . Minimum Masses from Linear Trends Alone 

In some cases, we find that a secular increase or de¬ 
crease in the observed radial velocities is present (a “lin¬ 
ear trend”), which is presumably a small portion of a 
Keplerian signal from a massive companion, typic ally an 
outer planet, or a secondary star or brown dwarf dCrepp 
et al. 1120121 |2013a|bl |2014[ |Montet et al.||2014[ |Knutson 
et al. 2014). 


When the trend shows no curvature and we have no 
AO imagery to put limits on the mass and angular sep¬ 
aration of companions, we usually say very little about 
the companion beyond a minimum mass (and a maxi¬ 
mum luminosity from the fact that its spectrum did not 
complicate the RV analysis). The scenario that gives the 
minimum mass to a companion generating a linear trend 

O 

of a given magnitude is one that has e ~ 0.5, w = 90 , 
which produces a sawtooth-like RV curve with a long, 
nearly linear component for ~ 80% of the orbit with a 
brief, high-acceler ation componen t during periastron for 
the other ^ 20% (Wright (2006), and see top panel of 
Figure [2 for an example; other panels show other patho¬ 
logical cases with radically different periods and semi¬ 
amplitudes that mimic the same trend). 

The minimum mass of a planetary companion detected 
only by its strongly detected constant acceleration 7 , 
is thus derived by s olving the mass functio n for the 
minimum mass (e.g., Wright & Howard 2009) assuming 
P ~ 1.25r (where t is the span of the observations), 
e ~ 0.5, and K ~ ry: 


M, 


minimum 


(0.0164Mj,p) - 


i^) 


4/3 


7 


ms ^yr ^ 


Mq 


2/3 


(1) 




Time since first observation (yr) 



Fig. 1.— Fifty synthetic RV measurements made over 8 yr by 
an unlucky observer of a hypothetical system with Gaussian errors 
of 3 m s“^. Top: RV curve of a planet with M sini = 1.97Mjup, 

P = 10 yr, R' = 30 m s“^, e = 0.5 and a; = 90 . The observer 
might conclude, incorrectly, that they were seeing the effects of a 
distant exoplanet with P ^ 8 yr and ^ 30 m s“^. Middle: 
RV curve of a planet with Msini = 6.6Mjup, P = 11 yr, = 

345 m s“^, e = 0.97 and lj = 20 . Although the period P and 
magnitude of the observed trend are about the same as that in 
the top panel, the true semi-amplitude of the orbit is much larger. 
The inset illustrates the complete RV curve, with the same units 
as the main figure. Bottom: RV curve for hypothetical planet with 
Msini = 11.7Mjup, P = 200 yr, A = 65 m s“^, e = 0.5 and 

o 

uj = 270 . Although the magnitude of the observed trend is the 
same as that in the top and middle panels, the period in this case is 
much longer, while K is only modestly larger. The inset illustrates 
the RV curve, with the same units as the main figure, over more 
than a complete orbit. The box in the inset illustrates t he span 
of the main panel. All three panels are reproduced from |Wright| 
l |200^ . 
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3. REFINED ORBITAL PARAMETERS FOR SEVEN 
PLANETARY SYSTEMS 

Our sample includes many known planetary systems 
of interest because of the presence of a linear trend in 
the residuals indicative of an additional companion; some 
with known trends with significant curvature; and some 
known to have outer companions with poorly-constrained 
parameters. The first six, HD 24040, HD 66428, HD 
74156, HD 183263, HD 187123, and HD 217107 are G 
stars; the seventh GJ 849, is an M dwarf. 

Table presents the time span of sets of observations, 
the numoer of points from each set, the number of new 
points, the mean uncertainty in velocities from each set, 
and the offsets between instruments. 

3.1. HD 24040 



Wright et al.| (2007 1 
ion to the star H 


lU 


reported a substellar compan- 
^040 with a wide range of possi¬ 


ble periods (10 yr < P < 100 yr) and minim um masses 
(5 < Msini < 20Mjup). Boisse et al. (2012), combining 
velocities from HIRES, SUHHiE, and ELUDIE, deter¬ 
mined an orbit of 3668^171 days (corresponding to 10 yr) 
and a minimum mass of 4.01±0.49Mjup for HD 24040 b. 

(20121 also found a linear trend of 3.85l|^'2g 


Boisse et al. 


m s 




, indic ative of a third body in the system. 
(]( 2012 ) also investigated potential long-term 
06 
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correlation between SOPHIE measurements and stellar 
activity indices but did not find such behavior. 

We present an updated fit with more recent Keck- 
HIRES velocities, seen in Eiguref^and Table We use 
HIRES data and published SOPIHE and ELODIE data, 
so in our fit we applied jitter instrument-by-instrument. 
With 107 velocities in total, 47 of which are from HIRES, 
13 from SOPHIE, and from 47 ELODIE (Boisse et al. 


2012 1 , we find for the best-fit one-planet Keplerian model 
an rms of 13.62 m s“^ and xi of 0.93. HD 24040 b orbits 
at a semimajor axis of 4.637 ± 0.067 AU, correspond¬ 
ing to a period of 9.5 yr, making it a good Jupiter ana¬ 
log in terms of its orbit (however, its minimum mass is 
4.10±0.12Mjup). The linear trend is 1.8±0.4 m s“^ yr“^ 
(lower than reported in |Boisse et al. 2012), a minimum 
mass of at least 1.44 Mjup accordiiig to Equation ([^. 
Our fit for HD 24040 6 , with a period of 3490 ± 25 days 
and minimum mass of 4.10 ± 0.12Mjup, is in good agree¬ 
ment with the solution from Boisse et al. (2012). 

3.2. HD 66428 


Butler et al. (2006) announced HD 66428 b, a planet 
with H = 1973 ± 31 d (5.4 yr), e = 0.465 ± 0.030, and 
M sin(i) = 2.82 ± 0.27Mjup. We update the orbital pa¬ 
rameters with a total of 55 velocities from HIRES (see 
Figure!^. The original fit used 29 velocities taken with 
HIRES from 2000 to 2006. Our new fit adds 26 new data 
points through late 2013. Capturing two complete orbits 
of HD 66428 b, the fit has an rms of 3.14 m s“^ where we 
assumed a jitter of 3 m s“^ and xt of 0.96. We determine 
a period of 2293.9 ± 6.4 days, or 6.3 yr. We determine 
a minimum mass of 3.19 5 ± 0.066M,i„ n, which is more 
massive than reported in 
Given our larger set o: 


Butler et al. (2006). 
radial velocities, it is under¬ 
standable that our s olution does n o t mat ch with the so¬ 
lution announced in Butler et al. (2006). The final fit 


OJ 

3 

■a 



Date (year) 

(b) Residuals 

Fig. 2.— Radial velocity and Keplerian fit for HD 24040 b. Solid 
lines represent the best-fit Keplerian orbits. The fit includes a 
linea r trend of 1.8 di 0.4 m s“^ 


|2(a)|Keck, SOPHIE, and ELO DIE RVs overplotted by best-fit one- 
pianet Keplerian model. |2(b)|Residuals of the RVs with the best-fit 
one-planet Keplerian model subtracted. 


m S 


yr " (corresponding to a minimum mass for the 


outer companion of at least 1.77 Mjup, by Equation ([I)). 

We run the fit with no jitter and no trend in order 
to see the significance of the detected trend. Eor that 
case, xi is 52.56, and the rms of the residuals is 7.46 
m s“^. To compare, we found an rms of 3.14 m s“^ 
and xl of 8.23 for seven free parameters (including the 
trend) and no jitter. Given the improvement in the fit 
with a trend included, the trend is significant. We also 
note that the eccentricity of the orbit is large: 0.442 ± 
0.016. The trend may indicate that the outer companion 
has influenced the orbit of the b component. Eurther 
monitoring will determine the nature of the source of the 
trend (i.e., whether it is due to a stellar or planetary 
companion). 

3.3. HD 14156 


Naef et al. (2004) described the HD 74156 two-planet 
system as a 1.86 ±0.03Mjup planet in a 51.64 ± 0.011 day 
period with a 6.17 ±0.23Mjup outer companion in a 5.5 
yr orbit. Mult iple authors have suspected a third planet 
in the system. Barnes & Raymond (2004) predicted one 
based on th e Pack ed Flanetary bystem hypothesis, and 


finds a previously unreported linear trend of —3.4 ± 0.2 


Bean et al. (2008) claimed the discovery of a compan- 
ion with H = 336 da ys as the thir d planet. Based on 
analysis of RV jitter, Baluev (2009) questioned the va- 
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2000 2002 2004 2006 2008 2010 2012 2014 
Date (year) 


(a) HD 66428 b 



Date (year) 

(b) Residuals 

Fig. 3.— Radial velocity and Keplerian fit for HD 66428 b, with 
a trend of 3.4 ni s“^ vr~^ incor porated. Solid lines represent the 
best-fit Keplerian orbits. I3( a)l Ke ck RVs overplotted by best-fit 
one-planet Keplerian model. P(b)| Residuals of the RVs with the 
best-fit one-planet Keplerian model subtracted. 


lidity of HD 74156 “d” as a fa lse detection due to annual 
systematic errors from HRS. Wittenmyer et al. (20091 
con cluded that the third p lanet was unlikely to be real, 
and Meschiari et al. (2011) updated the system with fur¬ 
ther oEserrailons'md'Teached the same conclusion. 

Here, we combine 226 velociti es from CORALI E and 
EL ODIE (44 and 51 o bservations Naef et al.|2004 ), HRS 
(82 |Bean et al.||2008|), and HIRES (52j (see Figure li]). 
We apply a two-planet Keplerian model. We added jit¬ 
ter instrument-by-instrument, and our fit has an rms of 
11.03 m s“^ and xl of 0.97. We have captured at least 
two orbits of HD 74156 c, making our orbital solution 
more robust than previously reported solutions. Table 
lists the orbital parameters. HD 74156 c is one of the 
more massive planets we have examined, with minimum 
mass 7.997± 0.095Mjup. Both planets have large orbital 
eccentricities (e = 0.64 and e = 0.38 for b and c respec¬ 
tively). 

In Eigure we plot the Lomb-Scar gle periodogram 
(jScargle |1982f|Horne & Baliunas||1986| of the residuals 
to our best two-planet ht. There is no indication of any 
power at the period of the purported d component, a 
result which is consistent with prior refutations of this 
signal (indeed, our analysis here uses much of the same 
data as previous work on the topic). Indeed, there is no 
hint of significant power at any period, indicating that 
there is no detectable third planetary companion in this 



2000 2002 2004 2006 2008 2010 2012 2014 
Date (year) 


(a) HD 74156 system 



2000 2002 2004 2006 2008 2010 2012 2014 
Date (year) 


(b) Residuals 




Fig. 4.— Radial velocity and Keplerian fits for the HD 74156 
syste m. Solid lines represent the best-fit Keplerian orbits. 
|4(a)| CORALIE, ELODIE, HRS , and Keck RVs overplotted by best- 
tit two-planet Keplerian model. |4(b)|Residuals of the RV s with the 
best-fit two-planet Keplerian model subtracted. |4(c)|and|4(d)| the 

1-n T _ r_ .. TTTN /I -1 r-/-* 7- _-1 __ _ _j : _ 1 * * * * 


r'nrvrec fr^r 7A1 h anH n 
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Fig. 5.— Periodogram of the residuals to our best two-planet, 
five-instrument fit to the RV data for HD 74156. There is no indi¬ 
cation of significant power at any period, or of any power at all at 
336 days, the period of the purported but disproven d component. 
We have computed the 99% false alarm probability in this figure by 
calculating the highest peak in each of 10,000 such pe riodograms 
calculate d for synthetic data sets of RV residuals (e.g., (Howard et| 
|al.|2009| |. We calculated each of the 10,000 synthetic sets by ran¬ 
domly assigning the actual residuals (drawn with replacement) to 
each of the times of observations of the actual observations. In 
99% of cases the tallest peak had power below 13.5. 

system. 


3.4. HD 183263 


First reported by Marcy et al. (20051, the HD 183263 
system showed a residual linear tre nd in addition to a 
3.7Mjup planet in a 634-day period. |Wright et al.|(|2007l 
attributed the new and signific ant curvature m tlie r esid- 
uals to an outer companion. Wright et al. (2009) fol¬ 


lowed up and constrained the minimum mass (3.57 ± 
0.55 Mjup) and period (8.4 ± 0.3 yr) for the outer com¬ 
panion, HD 183263 c, to which we report an updated set 
of parameters. 

With 66 velocities from HIRES, we implemented a fit 
with an rms of 3.68 m s“^ and an assumed jitter of 3.2 
m s“^. Figure ^ presents the RV curves for the sys¬ 
tem as well as the residuals. The orbit for HD 183263 
c appears to have finally closed, and it is significantly 
closer to circular (e = 0.051 ±0.010) and has a longer pe¬ 
riod than the solution from Wright et al. (2009), which 
found e = 0.239 ± 0.64 and P ^ 8.5 yr. We find for HD 
183263 c, that P = 4684 ± 71 days, or 9.1 yr; M sin* is 
6.90±0.12Mjup. While our best fit orbital solution does 
not match well with the previous orbital solution, our 
solution resides comfortably wit hin the stable portio n in 
the Pc-Mc sin ic space found by Wright et al. (2009 see 
their Figure 3). 


3.5. HD 187123 
Butler et al. (11998) discovered 


HD 187123 6, a 
0.52Mjup planet in a 3-day orbit. After many years of 


continued monitoring of this system, Wright et al. ([2007 ) 
announced a long-period outer companion with P TTD 
yr and a min i mum mass between 1.5 Mjup and 10 Mjup. 
Wright et al. (2009) presented a solution that constrained 
the mass and period of an outer companion to within 
20%, with P = 10.4±1.2 yr and M sin* = 2.0±0.3 Mjpp. 
Figur e sho ws an updated ht with HIRES data. jNaefj 


et al. (2004|) provide ELODIE velocities; however, since 
they have significantly worse precision and do not add 
temporal coverage, we do not use them here. The 108 



2002 2004 2006 2008 2010 2012 2014 
Date (year) 


(a) HD 183263 system 



Date (year) 
(b) Residuals 



Phase 

(c) HD 183263 b 



Phase 

(d) HD 183263 c 


Fig. 6.— Radial velocity and Keplerian fits for the HD 183263 
syste m. Solid lines represent the best-fit Keplerian orbits. 

6(a) Keck RVs overplotted by best-fit two-planet Keplerian model. 

Residuals of the R Vs with t he best-fit two-planet Keplerian 
model subtracted. |6(c)|and|6(d)| the RV curves for HD 183263 b 
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Keck observations still cover multiple orbits of the plan¬ 
ets; assuming a jitter of 2.23 m s“^, we find an rms of 
2.66 m s“^. From our fit, the period of HD 187123 c is 
9.1±0.13 yr and the minimum mass is 1.818±0.035Mjup. 
HD 187123 c appears to be a Jupiter analog, although 
its orbit is somewhat eccentric at e = 0.280 ± 0.022. 


3.6. HD 217107 


Fischer et al. (1999) presented HD 217107 6 as a 1.27 


M.iun planet in a 7.12-day period. A few years later, Fis¬ 


cher et al. (2001) identified a linear trend in the residuals. 


whi ch w^ likely caused by an outer companion. [Vogt et| 
^ (2005) reported the first orbit for HD 217107 c, mod¬ 


estly constrai n ed at P = 8.6 yr and M sinf = 2.1Mjup. 
Wright et al. (2009) constrained the orbit and mass of 
HD 217107 c to almost within 10%, with P ~ 11.7 yr 
and the minimum mass ^ 2.6Mjup. 

As with the case of HD 74156, we also have data taken 
by different teams from several instruments, we added 
jitter instrument-by-instrument. In our fit, we use veloc¬ 


ities from Keck (128 observati ons), Lick dWrigh t et ah 
121), and CORALIE (63 |Naef et al.|2001[ ) to find 
a fit an rms of 10.29 m s“^ (see Figure 


Because the outer planet has only barely (apparently) 
completed an orbit, its orbital parameters may be es¬ 
pecially uncertain (and are particularly sensitive to the 
assumption that there is not a third, longer-period planet 
contributing significantly to the velocities). To explore 
the robustness of our derived orbital period of the c com¬ 
ponent as a function of its minimum mass, we have con- 
structed a map in P-M sin i space (a variety of what 


Knutson et al. (^2014) call “Wright diagrams”; see Pa- 


tel et al. ( |2L)07 ), Wright e t al. (200 ^ and siin i lar ap- 


et al. 


proaches tak en in, e.g., Dumusque et al. (2011); Boisse 
( 2012| )). In this inap all orbital parameters have 


been optimized (i.e., they are at their maximum likeli¬ 
hood in a minimum sense) for each pair of P^ and 
Me sin ic in the map (except for the offsets among the 
four instruments, which are held constant at their over¬ 
all best-fit values). 

Figurej^shows the contour map, revealing that the 
orbital period and minimum mass for HD 217107 c are 
well constrained with P = 14.215)')q'q 4® yr and M sini = 

4.511o;q2-^Jup- These uncertainties are roughly consis¬ 
tent with the uncertainties determined via bootstrap¬ 
ping, which yields P = 14.215 ± 0.06 yr and M sinz = 
4.51 ± 0.07Mjup. This validates our choice of stellar jit¬ 
ter for this star, since the contours in the x^ maps are 
sensitive to the choice of jitter, while the bootstrapping 
uncertainties are almost completely independent of it. 
We report the bootstrapping uncertainties in Table 1^ 

To test the importance of our assumption that there 
are only two planets contributing detectable accelera¬ 
tions to the star, we repeated our bootstrapping analysis 
with a model that includes an additional, linear trend to 
the data. Though there is no statistical need to include 
such a trend in our model, giving our model the free¬ 
dom to include one could, in principle, affect the best-fit 
parameters for the outer planet. Indeed, though the pa¬ 
rameters of the b component do not change significantly 
in this model (as expected given its high frequency) we 
find a slightly different best-fit with such a model, with 
Pc, Tp^c, and Kc all changing by 2-4 cr, resulting in a min- 






Fig. 7.— Radial velocity and Keplerian fits for the HD 187123 
syste m. Solid lines represent the best-fit Keplerian orbits. 

7(a) Keck RVs overplotted by best-fit two-planet Keplerian model. 
TfE} 3,esiduals of the R Vs with the best-fit two-planet Keplerian 
model subtracted. |7(c)|and|7(d)| the RV curves for HD 187123 b 
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imum mass for the outer companion of 4.37Mjup. The 
uncertainties on the parameters of the c component in 
the with-trend model are larger by a factor of 2-4, com¬ 
fortably including most of the parameter estimates from 
the no-trend model. We conclude that our choice not to 
include a linear trend does not have a large effect on our 
conclusions or parameter estimations. 

3.7. GJ 849 
3.7.1. Orbital Fit 

Unlike the other stars in this work, GJ 849 is an M3.5 
dwarf. Various studies of this star’s com position have 
all found si milar, super-solar abundances: [Rojas-Ayala] 
et al. ( 2012) find [Fe/H ] = 0.3 1T0.17 (from K -band fea- 
tures); 


On ehag et al. 


band); and [I'errien et a 


(2012) find 0.35±0.10 (using J- 
( |2012 ) found 0.31±0.12 (using 


AT-band). 

GJ 849 hosts the first planet discovered orbiting an 
M-dwarf with a se mi-major axis greater than 0.21 AU. 
Butler et al. (2006) announced GJ 849 b, with P = 5.16 


yr and minimum mass 0.82Mjup. At the time, there was 
evidence of a linear trend of —4.75 m s~^ yr~^ , indicative 
of a second companion. Bonfils et al. (2013) also fitted 


the system with one planet and a iinear trend of —4.0 


,-i 


,.-1 


m s^" yr^", adding their HARPS data to the published 
HIRES velocities. 

Stellar magnetic activity had to be ruled out as the 
source of the trend. 


(|2012|) mon- 
for 


Gomes da Silva et al. 

itored several M-dwarfs from the HARPS program 
long-term magnetic activity. For GJ 849, they saw mild 
correlation in our velocities with the Na I index data. 
However, the a mplitude was not large enough. [Mon-j 


tet et al. (2014) provided the first orbital parameters for 
GJ 849 c, finding Msini = 0.70 ± 0.31Mjup, and 
yr, and found no correlation between stel- 


on 

p = i9.3;iy 


lar magnetic activity and the long-period signal of this 
outer companion. 

Our fit, using 35 velocities from HARPS (Bonfils et al. 


2013) and 82 velocities from HIRES spanning from 1997 


through early 2014, has further constrained the orbital 
parameters of the GJ 849 system. We incorporate a jitter 
of 3 m s“^, and our fit has an rms of 3.72 m s“^. 

GJ 849 6 is a 0.911Mjup planet in a 5.27 yr period with 
an orbital eccentricity of 0.038. GJ 849 c is a 0.944 ± 
0.07Mjup planet in a 15.1±0.66 yr period with an orbital 
eccentricity of 0.087 ± 0.06. 

GJ 849 c has the longest robustly measured orbital 
semimajor axis of any planet orbiting an M dwarf dis¬ 
covered to date. Indeed, it has one of the longest well- 
measured periods of exoplanets orbiting any kind of star. 
Exoplanets with similar period a nd period uncertainties 


in the litera ture include 55 Gnc d (Marcy et al .|2002 |Endl 

^ Ma 


et al.|201 2): HD 166724 b a nd HD 219077 b (Marmier et 
ar|2U13j); and HD 139315 (Howard et al. 2010|) — but 
tHese all orbit s tars with M > O.HIWq and the two from 
Marmier et al.j show significant eccentricity. The exo- 
planet with the longest robust ly measured orbital pe riod 


is /3 Pictoris (P = 20.5_i ^ yr Macintosh et al. 2014). 


We estimated the model parameters for GJ 849 m two 
additional ways to check for consistency and robustness. 
It is unclear whether the bootstrap resampling proce¬ 
dure provides an accurate estimate of GJ 849 c’s orbital 
parameters. In particular, the poor phase coverage be- 
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(a) HD 217107 system 




Phase 

(c) HD 217107 b 
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Fig. 8.— Radial velocity and Keplerian fits for the HD 217107 
syste m. Solid lines represent the best-fit Keplerian orbits. 

8(a) Keck RVs overplotted by best-fit two-planet Keplerian model. 

Residuals of the R Vs with the best-fit two-planet Keplerian 
model subtracted. |8(c)|and|8(d)| the RV curves for HD 217107 b 
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14.0 14.1 14.2 14.3 14.4 

Period (yr) 


Fig. 9.— Best-fit 100 x 100 map for fixed values of and 
Me sin ic for HD 217107 c. This confirms that the period and 
mass are well-constrained. We have illustrated the contours of 
the la, 2a, and Sa (defined by = Xmm {2.30,6.17,11.8}) 
confidence le vels, based on for t he number or degrees of freedom in 
the problem (|Press et aL|2002||. The center and la limits in both 
parameters are consistent witti the bootstrapping uncertainties for 
these parameters. 


fore 2001 (see Figure 10) results in several clear outlier 
models in the joint parameter distributions. 

In our first check for accuracy in the parameters and 
uncertainties, as with HD 217107 c, we constructed a P- 
M sin i map to confirm that the orbital period of GJ 
849 c is well constrained (assuming no additional plan¬ 
ets and a stellar jitter of 3 m s“^), despite having just 
completed an orbit, and find that the 68% confidence in¬ 
terval contours corresponds to uncertainties in P of less 


sponi 

than 5%. As Figure |ll] shows, the map uncertain¬ 
ties in minimum mass are 0.07 Mjup, exactly consistent 
with our bootstrapping errors; the uncertainties in pe¬ 
riod are ~ 1.1 yr, which is larger than the bootstrapping 
errors of 0.66 yr, probably because the contours are 
asymmetric. 

As a second check, we turn to a Bayesian approach 
for performing parameter estimation via Markov chain 
Monte Carlo. We adopt the usual broad priors for Ke- 
plerian orbital parameters and likelihood assuming un¬ 
correlated, Gaussian measurement errors with dispersion 
based on the quadrature sum of the reported measure¬ 


ment uncertainties and an unknown jitter term (Ford 


20061. Given the potential for mutual planetary inter¬ 
actions, we apply RUN DMCj^ a well tested code that 
combines n-body integration with differential evo lution 
Markov chain Monte Carlo (jNelson et al. 2014). Al¬ 
though the GJ 849 planets are well approximated by 
Keplerian orbits, the differential evolution proposal in 
RUN DMC is much more efficient than a traditional ran¬ 
dom walk MCMC for dealing with correlated parame¬ 
ters, which are often present in the parameters for long- 
period companions, and so by using RUN DMC we do 
not have to fine tune a proposal distribution. We find 
that the marginal posterior probability distribution for 
Pc has 68% of its mass within 0.74 yr of the median pe¬ 
riod of 15.1 yr, only slightly larger than the uncertainty 
estimated from the bootstrap. 


^ We used t he Keplerian parameter priors given in |Nel-| 
|son et al.| (|2014|l, and the algorithmic parameters nchains=liUO, 
ngen=ruu,0uu, (7l,=0.01, and MassScaleFactor=1.0. 


The similarity of the parameter uncertainties from all 
three methods verifies that the orbit of GJ 849 c is well 
constrained and validates the BOOTTRAN and map 
approaches (in this case) and our choice of jitter. We use 
the more conservative contours do determine param¬ 
eter uncertainties in Table |3l 

3.7.2. Stability 

Because th is system is not “highly hierarchical” 
(Wright 2010) in mass or orbital period, we have per- 
lormed n-body simulations to establish the dynamical 
stability of our orbital solutions. The 1000 BOOTTRAN re¬ 
alizations of the GJ 849 RV data are used to determine 
parameter uncertainties is associated with a complete set 
of Keplerian orbital parameters for the two planets (P, 
e, w, AT, and Tp for each planet, plus an overall RV offset 
7 and two offsets among the three RV data sets). All 
of these realizations returned reasonable fits, indicating 
that the fitting procedure did not fail in any case. 

We performed long-term dynamical integrations for 
all 1000 fits to these realizatio ns of the data us ing the 
MERCURY symplectic integrator (Chambers 1999). Each 
simulation runs for 10^ orbits of the inner-most planet 
(^ 2 X 10^° days). This integration timescale is short 
relative to the lifetime of the star but sufficiently long 
enough to show a significant fraction of our models un¬ 
dergo an instability, described below. 

An instability occurs if at any point during the inte¬ 
gration either planet crosses the other’s Hill sphere or 
either of the planets’ semi-major axes change by more 
than 50% of their initial value. 

None of our models resulted in a collision over the 
course of the integration. However, we find 67 models 
undergo the second listed mode of instability (| [ofinai — 
ainitiai]/ainitiai |> 0.5) when the periastron passage of GJ 
849 c is less than 3.5 AU (Figure 12). The instability 
times are logarithmically uniform from ~10 to ~ 10^ yr. 

We removed the unstable BOOTTRAN realizations from 
our calculations of the uncertainties in the orbital pa¬ 
rameters we report in Table 

4. A 2Mjup PLANET AROUND HD 145934 

We here announce a new long-period planet orbit- 


(Takeda et al. 

2007 

). This star was not known to be 

a giant when t 

he Caiitornia Planet Survey began mon- 


itoring it in 1997 at Keck Observatory. Since then, its 
log(g) value and mass from Takeda et al. (2007) indicate 
that it is a giant. Visual inspection of the gravity sensi¬ 
tive sodium and magnesium lines confirm this diagnosis. 
Radial velocities of HD 145934 show a clear sinusoidal 
modulation of planetary amplitude upon a large linear 
trend, indicative of a stellar binary companion. 

In our analysis of the 75 HIRES velocities for HD 
14593 4, we n ote the slight overall curvature present (see 
Figure|13(d)|). To account for the curvature using RVLIN, 
which (at the moment) only accommodates purely linear 
trends, we treated HD 145934 as a two-companion sys¬ 
tem, with the outer companion having a very long (60 
yr) orbital period and circular orbit. There is not enough 
information in our time series for the resulting orbital pa¬ 
rameters of the outer companion to be meaningful, but 
this approach provides us sufficient flexibility to fit out 
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the low-frequency power contributed by the binary com¬ 
panion. Equation Q constrains the minimum mass of 
the companion to be at least 21Mjup. 

To determine the effects of modeling the ostensible stel¬ 
lar companion with our choice of orbital parameters on 
the planet’s parameters, we checked first the impact of 
letting eccentricity be a free parameter. The best-fit ec¬ 
centricity is close to circular ( 0.05), so our choice of 
fixed e = 0 is not strongly affecting our analysis. We 
also changed the (fixed) period of the stellar companion 
to take values between 50 and 80 yr (guesses outside the 
range returned poor fits, but given the nonlinear nature 
of the problem this does not necessarily reflect an actual 
upper limit to the companion’s period). We found that 
the choice of period did not have significant impact on 
the parameters of the planet. For example, the best-fit 
values for the period of HD 145934 b varied on the order 
of 10 days for different outer companion periods. The 
minimum mass varied on the order of 0.1 Mj^p. These 
differences are all well within Icr of our presented set of 
parameters. We conclude that our modeling of the outer 
companion is sufficiently flexible to have no important 
effects on our estimates of the planet’s orbital parame¬ 
ters. 

Given that the rms of the residuals to the fit without 
stellar jitter is 7.83 m s~^. we assume a st ellar jitter of 
7.5 m s“^ in our fit. Hekker et al. (2006) performed a 
survey of stable K giants with jitters lower than 20 m 
s“^. The most stable of that sample range between 6 
and 15 m s“^, so our choice of jitter is reasonable and 
also consistent with the residuals. The residuals to the 
resulting best-fit Keplerian model have rms of 7.80 m s“^ 
and xl of 1-05. We find that HD 145934 b has a period 
of 7.48 ± 0.27 yr, an orbital eccentricity of 0.053^Q gg 3 , 
and a semi-amplitude of 22.9±2.6 m s“^. The minimum 
mass of the planetary companion is 2.28 ± 0.26 Mjup. 

The presence of curvature in the binary companion’s 
orbit implies that either it is highly eccentric and near 
periapse, or that we have observed a nonnegligible por¬ 
tion of its orbit. The latter is more likely, and implies 
that its orbital period is a several or dozens of decades, 
not millennia. 


5. DISCUSSION 

Our analysis of 13 exoplanets uses recent Keck-HIRES 
radial velocities and other published data. We see that 
there is need for follow-up work, as in the cases of GJ 
849 and HD 145934 for better constraints and further 
analysis. In the instance of HD 66428, whose residuals 
show a previously unseen linear trend, we will monitor 
for the completion of orbits or to rule the companion out 
as a planet. 

We have reduced the uncertainties in the parameters 
for many planets. The up-to-date HIRES data allowed us 
to place upper limits or constrain several orbits. From 
our sample, we identify two planets as Jupiter analogs 
around Sun-like due to similarities in semimajor axis (5.2 
AU): HD 24040 b and HD 187123 c, although both are 
much more massive than Jupiter, and the latter’s orbit 
is somewhat eccentric. We have discovered a new planet, 
HD 145934 b, and its host star’s residuals show curvature 
whose velocity semiamplitude is indicative of a probable 
stellar or brown dwarf companion. 


We confirm GJ 849 c, and find that it is the planet 
with the longest known period around an M dwarf so 
far. GJ 849 is a rare system in that it is a multi-giant- 
planet system around an M-dwarf. In all of our multi¬ 
planet systems, the inner planet is less massive, though 
this fact is certainly influenced by the soft decrease in 
semiamplitude with orbital distance (K oc HD 

66428 may be a case where the planet’s high eccentricity 
and the presence of a linear trend in the system are signs 
th e outer companion has affected the inner planet’s orbit. 


Kane et al. (2014) found in the case of HD 4203. 


as 

All of these systems, but perhaps e specially the “highly 
hierarchical” systems (|Wright|2010|) HD 187123 and HD 
217107, will be valuable tor reconciling observations and 
the theory of planetary migration. These two systems 
are at present the only known examples of systems con¬ 
taining a hot Jupiter (gas giant with P < 10 days and 
Msini > O.lMjup) and a very-long-period planet (P > 5 
yr) with a well determined orbit. In both cases, the outer 
planet is ~ 3 times the mass of the inner planet, and there 
is no evidence of other planets in the system. 

There are only two other systems with hot Jupiters and 


iButler et al. 

2006 Wright et al. 

2009 

) and 

Bakos et al. 

2009 

Winn et al. 2 

OT^ 

'Lhe 


HAT-P-13 

former case remains anomalous in that the innermost 
planet is the most massive, with M sin* = 3.9Mjup (the 
outermost planet has Msini = 0.6Mjup and P = 2.6 yr; 
there is also a third, intermediate planet in the system). 
The latter case has an especially high mass ratio, having 
a highly eccentric Msini > 14Mjup outer planet and 
an inner, transiting planet with M = O.SdM .Tim. We 
know from bot h RV studies (Wright et al.||2009) and the 
Kepler results (Latham et aF 21)11) that “hot Jupiters 
are lonely”, at least when it comes to companions within 
^ 1 AU. Gontinued long-term monitoring of other hot 
Jupiters will establish whether they have freque ntly have 


cold friend s” at larger orbital distances (e.g., Knutson 

eTaTjlMit . 
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Fig. 10.— Radial velocity and Keplerian fits for the GJ 849 
system . Solid lines represent the best-fit Keplerian orbits. 


10(a)| Keck RVs overplotted by best-fit two-planet Keplerian model. 

Residuals of the RV s w ith the best-fit two-planet Keplerian 
model subtracted. |10(c)| and |10(d)| the RV curves for GJ 849 b 
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Fig. 11.— Contours at Itr, 2 (t, and Scr (defined by = Xmin + 
{2.30,6.17,11.8}) confidence levels, as appropriate for the number 
of degrees of freedom in the problem | |Press et al.|[^002| | for the 
orbital parameters of GJ 849c. The period and mass of the c 
component appear to be well constrained it better than 5%. 





5.0 



0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 
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Fig. 12.— Best fit minimum mass and periastron distance for 
all of the 1000 BOOTTRAN realizations of the GJ 849 system (which 
were generated through a bootstrapping procedure performed on 
the radial velocity data). Red x’s represent fits which proved to 
be unstable in fewer than 10^ orbits. There is a clear boundary 
at periastron distance ~ 3.6 (AU) between stable and unstable 
configurations. We have not used the unstable realizations in de¬ 
termining parameter uncertainties for the GJ 849 system. 




Fig. 13.— Radial velocity and Keplerian fit for HD 145934. Solid 
lines r epresent the best-fit Keplerian orbit. 

|13(a)|Keck RVs overplotted by best-fit two-planet Keplerian model 
tnat i ncludes a long-period companion to account for curvature. 
|13(b) Residuals of the RV s with th e best-fit two-planet Keplerian 
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TABLE 1 

Summary of Radial Velocity Data (Number of Observations) and Mean Uncertainties 



Instrument 

Mean Unc. (m s ^) A^obs Span (yr) 

A^new 

Offset from Instrument 1 (m s ^) 

HD 24040* 

1 

SOPHIE 

4.27 

13 

2008-2010 



2 

ELODIE 

11.48 

47 

1997-2005 


-53.19 ±4.37 

3 

HIRES (pre-upgrade) 

1.35 

20 

1998-2004 


-44.73 ± 3.80 

4 

HIRES (post-upgrade) 

1.56 

27 

2004-late 2013 

22 

-34.28 ± 1.69 




HD 66428. 

Jitter = 3 m s~^ 



1 

HIRES (pre-upgrade) 

1.22 

22 

2000-2004 



2 

HIRES (post-upgrade) 

1.05 

33 

2004-late 2013 

26 

2.8 ± 1.8 

HD 74156* 

1 

CORALIE 

8.52 

44 

2001-2003 



2 

HRS 

8.34 

82 

2004-2007 


-32.5 ± 3.07 

3 

ELODIE 

12.74 

51 

1998-2003 


32.16 ±2.82 

4 

HIRES (pre-upgrade) 

1.99 

9 

2001-2004 


47.02 ± 2.6 

5 

HIRES (post-upgrade) 

2.87 

43 

2004-late 2013 

31 

66.73 ±2.6 




HD 145934. 

Jitter = 7.5 m s”-*^ 



1 

HIRES (pre-upgrade) 

1.22 

44 

1997-2004 



2 

HIRES (post-upgrade) 

1.00 

31 

2004-early 2014 

75 

14.99 ± 9.70 




HD 183263. 

Jitter = 3.2 m 



1 

HIRES (pre-upgrade) 

1.6 

31 

2001-2004 



2 

HIRES (post-upgrade) 

1.23 

11 

2004-mid-2013 

24 

4.04 ±4.4 




HD 187123. . 

Jitter = 2.23 m s“^ 



1 

HIRES (pre-upgrade) 

1.22 

64 

1997-2004 



2 

HIRES (post-upgrade) 

1.19 

46 

2004-mid-2013 

40 

-1.58 ±0.88 

HD 217107* 

1 

Hamilton 

4.709 

121 

1998-2007 



2 

CORALIE 

9.175 

63 

1998-1999 


-402.3 ± 1.74 

3 

HIRES (pre-upgrade) 

1.414 

63 

1998-2004 


17.6 ± 1.38 

4 

HIRES (post-upgrade) 

1.022 

68 

2004-late 2013 

31 

25.68 ± 1.87 

GJ 849. Jitter = 3 m s ^ 

1 

HIRES (pre-upgrade) 

3.19 

24 

1997-2004 



2 

HIRES (post-upgrade) 

1.48 

58 

2004-late 2013 

3 

-4.24 ± 2.24 

3 

HARPS 

1.05 

35 

2003-2008 


18.36 ±2.1 


Note. — For each star, we list the chosen jitter, instruments, the mean uncertainty corresponding to the data from each 
instrument before jitter is applied, number of observations, date range of observations, the number of new points this work has 
added, and offsets from instrument 1. In a few cases, marked by an asterisk (*), we add jitter instrument-by-instrument. We split 
data from HIRES by the 2004 upgrade. 

For references of data, see text. 



TABLE 2 

List of Stellar Parameters 


Updates to Long-Period Planets 


0 ? 
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_ __ X o 
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CTJ u 




O iO 00 t> CO o 
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) ^ CO t> X 

I CO CO X CM CO o 

; O CM ^ ^ O --H S 

I . . K, . . . CCl 

^ ^ ^ ^ ^ ^ 


CO O I-H 
O ^ X 
CM X I-H 
CO CO d 


1 4.361(60) 

1 4.490(60) 

1 4.259(60) 

3.23(6) 

1 4.403(60) 

1 4.359(60) 

1 4.541(60) 

^ ; 

^ ^ ^ 

^ : 
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parentheses, has the same place value as that of the quantity. For example, “0.460(30)” is equivalent of “0.046 ± 0.030”, and “5898(44)” is equivale nt of “5898 zfc 
44”. Unless stated otherwise: the values in Columns 2-6 are from various sources collated in the SIMBAD Astronomical Database; distances are from|van Leeuwen 
<20071; dat a in Columns 8—11 are from [Valenti 8z Fischer] j2005] ; data in Column 12 are from |Takeda et al.| | f2007] ; and we calculated the stellar radii based on 
[Torres et al.|(201^ . 

^ For GJ 849, T^ff is from [Rojas- Ayala et ai~|j201^ , [Fe/H] is from [^rrien et al.|(^012] , while M* is from [Montet et al.|(^014] . 
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TABLE 4 

Radial Velocities Measured for HD 
24040 


BJD 

-2440000 

RV 

(m s~^) 

itlfJ 
(m s“^) 

Tel 

10838.77321 

-95.2 

1.5 

HIRES 

11043.11965 

-85.4 

1.5 

HIRES 

11072.03904 

-87.0 

1.5 

HIRES 

11073.00232 

-87.8 

1.1 

HIRES 


Note. — This table is available in its entirety 
in machine-readable form. 



